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Abstract. 

The effects of a magnetic field (h) and a space modu- 
lation (5) on the magnetic properties of a one dimensional 
antiferromagnetic-ferromagnetic Heisenberg spin-1/2 model 
have been studied by means of numerical exact diagonal- 
ization of finite size systems, nonlinear sigma model and 
bosonization approach. The space modulation is considered 
on the antiferromagnetic couplings. At 8 = 0, the model 
is mapped to a gapless Liittinger liquid phase by increasing 
the magnetic field. However, the space modulation induces 
a new gap in the spectrum of the system and the system ex- 
periences different quantum phases which are separated by 
four critical fields. By opening the new gap a magnetization 
plateau appears at ^M sat . The effects of the space modu- 
lation are reflected in the emergence of a plateau in other 
physical functions such as the F-dimer and the bond dimer 
order parameters, and the pair-wise entanglement. 
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1. Introduction 

The effects of a magnetic field on the magnetic properties of low-dimensional quantum 
magnets at zero temperature has attracted much attentions in recent years. One 
of the very interesting phenomena is the appearance of a magnetization plateau in 
the magnetization curve. This behavior can be viewed as an essentially macroscopic 
quantum phenomena, and has gained much attentions, recently. When a plateau is 
appeared the energy gap is opened, which can be in some senses regarded as a kind of 
generation of the Haldane conjecture pp. In a seminal work, M. Oshikawa, M. Yamanaka 
and I. Affleck studied the magnetization of a general class of the Heisenberg spin chains 
at the presence of a magnetic field. They showed that the plateaus can be appeared 
when the magnetization per site m is topologically quantized by n(s — m) = integer, 
where s is the magnitude of the spin, and n is the period of the ground state determined 
by the explicit spatial structure of Hamiltonian [2]. 

The bond alternating Heisenberg spin- 1/2 chains which are obtained by a space 
modulation in the exchange couplings [31 HI El El [7J El |9] are a particular class 
of the low- dimensional quantum magnets to observe the magnetization plateau at 
zero temperature. The bond alternating Antiferromagnetic- Ferromagnetic spin-1/2 
chains have a gap in the spin excitation spectrum and reveal extremely rich quantum 
behaviors in the presence of an external magnetic field pHl El EEl [13]. By turning 
the magnetic field, the excitation gap is reduced and reach to zero at the first critical 
field. Simultaneously, the magnetization start to increase up to its saturation value, 
0.5 at the second critical field. More enhancement of the field re-opens the gap and 
the saturation plateau appears in the magnetization curve. These models have only 
two plateaus at zero and saturation values of the magnetization. It has been found 
that a space modulation in the exchange couplings can affect on the behaviors of 
the field induced magnetization. For example the bond alternating Ferromagnetic- 
Ferromagnetic- Antiferromagnetic (F-F-AF) trimerized Heisenberg spin-1/2 chains have 
exotic behaviors by changing the magnetic field [2]. This model can be realized in 
the Cw-compounds such as 3CuCl2-2dx. The magnetization has a plateau at = |, 
where M and M sat are the magnetization and its saturated values, respectively |T3l [TU] . 
The mid-plateaus have also been appeared in the bond alternating Ferromagnetic- 
Antiferromagnetic- Antiferromagnetic (F-AF-AF) trimerized Heisenebrg spin-1/2 chains 
[T7j . It has been also shown that the static structure factor dose not vary with the 
external magnetic field at the plateau state [T8] . 

The other examples are the tetrameric bond-alternating Ferromagnetic- 
Ferromagnetic-Antiferromagnetic-Antiferromagnetic (F-F-AF-AF) Heisenberg spin 
models. A realization for the spin-1/2 chain with F-F-AF-AF alternations is the com- 
pound Cu(3 — Clpy)2(N 3 )2(3 — Clpy = 3 — C hloroyridine) [19]. In such a spin-| system, 
there is a gap from the singlet ground state to the triplet excited states in the absence 
of a magnetic field. At — h a plateau is appeared in the magnetization curve 
by applying a magnetic field [20]. The temperature dependence of the magnetization, 
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magnetic susceptibility and specific heat have been also studied in the mid-plateau state 
of this model by means of transfer-matrix renormalization group method [21] . 

In this paper, by considering a different bond alternating Heisneberg spin- 1/2 chain, 
we study the other physical properties of the mid-plateau state such as the string 
order and pair-wise entanglement, F-dimer and bond-dimer order parameters. The 
Hamiltonian of the model is 

N/2 

H = — Jf • S2J+1 

N/2 N 

+ Jaf + (-1)^-1 • S 2j -hJ^Sj, (1) 

i=i 3=1 

where Sj is the spin-| operator on the j-th site. Jp and J^ F = Jaf(1 ± 8) denote the 
ferromagnetic (Jp > 0) and antiferromagnetic (Jaf > 0) couplings, respectively. 5 is 
the space modulation parameter and h is proportional to the external magnetic field. 
The unit cell of the model has been shown in Fig. (|lp. The bond alternating spin- 1/2 

— -9-—J- — » — j- F - — m -•- — 

Figure 1. (Color online.) Schematic representation of a tetrameric spin chain 

AF-F chains are good candidates for studying the Liittinger liquid phase. A suitable 
realization of the bond alternating AF-F spin chain is the Cu based (CH 3 )2NH 2 CuCl3 
compound. Linked-cluster calculations and bulk measurements show that DMACUCI3 is 
also a realization of the spin-1/2 alternating AF-F chain with nearly the same strength 
of antiferromagnetic and ferromagnetic couplings [22J. A space modulation changes the 
AF couplings alternatively and influences on the quantum properties of the system. 

Using the exact diagonalization of a finite size system and employing the conformal 
field theory it has been already found that for 5 7^ a magnetization plateau exist at 
half of the saturation value [23]. Moreover, the width of the mid-plateau is dependent 
to the values of S and the critical exponent for the plateau width is obtained by the level 
spectroscopy method. In the present work the different properties of the model such as 
F-dimer and bond-dimer order parameters, string order parameter, and the pair wise 
entanglement of the plateau state are studied. 

The outline of the paper is as follows. In the forthcoming section we discuss the 
model in the strong AF coupling limit and derive an effective spin chain Hamiltonian. In 
section III we summarize the results of the analytical field theory studies. In section IV, 
the results of exact diagonalization Lanczos method are presented. Finally, we discuss 
and summarize our results in section V. 
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2. Effective Hamiltonian 

The easiest way to obtain the effective model, which manifestly displays mechanism 
for generation of the sequence of new scales is to start from the limit of noninteracting 
block of pairs Jp = and strong magnetic field h ~ Jaf [21] • In the limiting case of 
the strong AF coupling Jaf 3> Jf an d Jaf 3> 5 Jaf the model can be mapped to an 
effective spin chain Hamiltonian [12]. At Jaf ^ Jf, the system behaves as a nearly 
independent block of pairs. Indeed an individual block of spin pairs may be in a singlet 
or a triplet state with the corresponding eigenvalues given by 

E ± (S) = ~^Jaf, 

E± {T l ) = -j\ F - h, E ± (T ) = -Jaf, 



E ± {T_ 1 ) = \j± F + h. (2) 

When h is small, the ground state consists of a product of pair singlets. As the magnetic 
field h increases the energy of the triplet state |Ti) decreases and at h — J^p forms 
together with the singlet state, a doublet of almost degenerate low energy states, split 
from the remaining high energy two triplet states. Thus, for a strong enough magnetic 
field we have a situation when the singlet \S) and triplet |T\) states create a new effective 
spin r = 1/2 systems. On the new singlet-triplet subspace and up to a constant, we 
easily obtain the effective Hamiltonian 

j N/2 

2 ^ j j+1 j j+1 2 j 

N/2 N/2 
0=1 3=1 

where /ig = h — J A f + ^ and = 5 Jaf- Thus, the effective Hamiltonian is 
nothing but the XXZ Heisenberg chain, with anisotropy A = 1/2 in a uniform and 
staggered longitudinal magnetic field. The full phase diagram of this model has been 
investigated by F. C. Alcaraz and A. L. Malvezzi [25] and the nature of the ground 
state phase transition has been pointed in Ref. |26j. The other properties also could 
be found in recent works such as Refs. [2J [271 ESI EH] • To find clearer picture for the 
low energy spectrum of the effective model ([3]), using the numerical Lanczos method we 
have computed the energy gap of the effective model. The numerical results, have been 
shown in Figj2]for a chain of lengths iV = 16, 20 and the exchange parameters Jf = 1.0, 
Jaf — 9/2, and 5 = 1/9. As it is seen from Figj2, at h e J^ = or h = Jaf — if, the 
effective model reduces to the XXZ spin chain in a staggered longitudinal magnetic field. 
The ground state of this model is in the gapped Neel phase. By changing from 
zero i.e decreasing (increasing) h from the value h = Jaf ~ if; the energy gap decreases 
and reaches to zero at finite values of the magnetic field h e J^~ {h e J^ ,+ ). Consequently, 
there is a plateau in the magnetization curve of the effective model at < h e J^ < h e J^ + . 
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By more decreasing (increasing) of h, the energy gap remains zero dawn(up) to critical 
field {h e J^' 2 ) and, after these gapless regions it is re-opened and behaves linearly 

with the magnetic field Hq . Correspondingly, in these regions the magnetization of the 
effective model starts to change from zero and is saturated at the critical field (/i^' 2 ). 

3. Field theory predictions 

In this section we will study the model ([T]) in the language of two continuum field 
theories. We have employed the nonlinear a model and Bosonization approaches to 
obtain the energy gap, the critical fields and magnetization plateaus of the model. 

3.1. Nonlinear a model 

The 0(3) nonlinear a model (NLaM) is a semiclassical approach which is based on 
the certain properties of a field theory in (1+1) dimension. In 1983 Haldane predicted 
theoretically the existence of a finite gap between the ground state and the first excited 
state of antiferromagnetic Heisenberg integer spin chains [T]. He also conjectured that 
the half-odd integer spin chains are gapless. The Hamiltonian of a homogeneous spin 
chain is mapped to a 0(3) NLctM with an additional topological term. The topological 
term is if spins of the chain are integer and tt if they are half-odd integer. 

Following we will investigate the low temperature behaviors of the alternating spin 
chain (jl|. For convenient on the incoming calculations let us write the Hamiltonian (jl]) 
as the following form 

N/4 
3=1 
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N/4 

+ Jaf £[(1 + 5)Sg ) _ 1 ■ Sg> + (1 - <$)sg> • S« J 

j 

(4) 

where < 5 < 1. Using the spin coherent states representation for the spin operators, 
i.e S = sSl we can write the Hamiltonian in terms of the classical spin vectors. By 
introducing three classical fields n, L, and A the classical vectors are written as the 
following forms: 

(i) - n fl | fl T l 2 W 2 _L_ a T 

- (-1)^A 2J . 1; 

ng =n 2 ,(l-|%| 2 )V 2 + ^ + (-l)^A 2j , 

J s s s 

(5) 

where i = 1,2 and a is the lattice constant. Typically, for a quantum antiferromagnetic 
Heisenberg chain it is convenient to write the spin vectors in terms of a unimodular 
Neel field (n) and one ferromagnetic canting field (L). However for an alternating AF-F 
Heisenberg spin chain, the main difficulty arises from the fact that it is inhomogeneous 
and writing a continuous action out of the discrete Hamiltonian is not a trivial task 
[30] . Using another field which describes the variation of the Neel fields inside the two 
spin-blocks such as A we can map the Hamiltonian Q to the following 0(3) NLcrM: 

cL pLn 



1 



22 



o jo 



A=— / dxdx [(d x n) 2 + (d X0 n) 2 ] - iOW : (6) 



where W is the winding number, Lt = c/3, xq = cr, and we have considered the case 
h = 0. By defining a = j^, the coupling constant g, and the velocity of spin excitations 
c are given in terms of a and S as: 



2(l+a) 

9 ~ 

I _ §1 _ 



4 2(l+a) 



1/2' 



6 2 1 X V2 



^ = ^a,(l- T -— ^) . (7) 



The topological term is obtained as follow: 

s5i 
■ 1 + a 



e = 2,( 7T ^L ? ), (8) 

2 

where . For a = or 5 = the topological term is zero and the model is always 
gapped and the spin excitations velocity is ^Jafo-s. Moreover, the a = 1 and 5 = is 
the special case of AF-F Heisenberg spin chains which the topological term is always 0. 
This result is in good agreement with the result presented in Ref. [30J. For any other 
values of a and 5, the Hamiltonian is mapped to a non-integrable 0(3) NLaM with 



Space modulation on ID alternating spin- 1/2 model 



7 



in the interval [0, 7r]. This model is also gapped and the gap value is dependent on a 
and S. 

Now, let us study the model in presence of the magnetic field. Since, the model 
is always gapped in the absence of a magnetic field we are allowed to consider 6 = 0. 
The Hamiltonian Q is mapped to the following NLcrM which the effect of the magnetic 
field is clearly seen: 

A =jgj d 2 x[(d x n) 2 + (d Xo n -*-hx n) 2 ], (9) 

To get more physical insight from the effects of the magnetic field, it may 
be more physically transparent to work with two fields (9, 0) where n = 
(sin 9 sin 0, sin 9 cos 0, cos 9). Here 9 is co-latitude and is azimuthal angles. Selecting 
magnetic field in z direction, the action ^ takes the following form; 

A=^- f d 2 x{(d^9) 2 + sm 2 9 f ' * 1 
^9 J 



'■ 2 + J_(0--^ 



cf 



}, (10) 



A = — I I d 2 fi(-h 2 - 0<9j0 - 2ihd xo c 



where 9 is the angle of h and n. In the term, h x n, the magnetic field induces a 
hard-axis anisotropy. In other word, the magnetic field try to align all spins with n in 
the plane normal to h. Thus for high field regime, h > |0'|, the deviation of Neel field n 
from the plane is small. Thus an expansion to quadratic order in •& is valid. The action 
is 

1 

+ («V) 2 - $ 2 [(^ o 0) 2 -h 2 - 2ihd x j] ) , (11) 

where h = h/c and the fluctuations have been separated to the in-plane and out-of-plane 
fluctuations. 

Using a spin stiffness analysis, 1/N expansion and a renormalization group 
approach, it has been already computed the magnetization and spin correlation functions 
of a spin ladder in an applied magnetic field [3TJ [32]. The magnetization (£\ is given 
by M = jf^, where F = —jjlnZ is the Helmholtz free energy and Z = J V-dVcfye'^. 
The separable nature of fluctuations allows us to give the results of M = M Q + Mj 
which is a summation of both out of plane and in plane contributions. At low enough 
temperatures and small value of c/h one can find that the out of plane contribution is a 
constant, M Q = 1/2, which is correspond to a uniform state. The in plane contribution 
has two terms, one is linear in the field and the other has a sawtooth form (See Ref. 
|32j). Totally, Mi is a step-like form which the width of the steps scales as 1/N. Thus 
the sawtooth form is the finite size corrections and, in the thermodynamic limit iV — > oo 
the in plane magnetization is only linear. 

For our tetramerized spin chain the magnetization is given by 

M--1 + — . (12) 

Application of an enough high magnetic field will cause the spin alignment, or saturation 
with a maximum magnetization, M s = s. This effect should be considered in the NLcrM 
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approach by a lagrange multiplier. At zero magnetic field the system is gapped and 
the magnetization is zero. Tuning magnetic field decreases the gap till first critical 
field, h Cl where magnetization start to increase. At the second critical field, where the 
magnetization saturated, the gap reopens and spins are fully aligned in the magnetic 
field direction. 

The critical fields which are attained by means of NLcxM are as follows: 

h ci = gc, 

3 

h C2 = -gc. (13) 
3.2. Bosonizaion 

In this section we concentrate our attentions to the low-energy and long wavelength 
excitations by using bosonization language. Let us consider the Hamiltonian ([3]). In our 
analysis of the model ^ we closely follow the route developed in the Ref. [27]. In the 
absence of both magnetic fields, Hq = and = we have the XXZ spin-1/2 chain 
with ferromagnetic coupling and anisotropy parameter A = 1/2. This model is critical 
and the long-wavelength excitations are described by the following Gaussian theory; 



n = -j dx[{d x ey + (d x( f>)% (i4) 

where 8(x) and <j){x) are dual bosonic fields, d t <p = vd x 9, and satisfying the following 
commutation relations 

[<f>(x),9(y)]=ie(y-x), (15) 

[<Kx),6(x)] = ± (16) 

v is the spin excitation's velocity and fixed by the Bethe ansatz solution as 

where Liittinger parameter K (inverse of Bosons radius [33]) is given as a function of A 

as 

K= — t -. (18) 

2 arccos A 

Thus K is increased monotonically along the XXZ critical line — 1 < A < 1 from its 
minimal value K — 1/2 (for the isotropic antiferromagnetic A = — 1 case) to unity 
at A = (for the noninteracting case) and goas to infinity at A = 1 which is the 
ferromagnetic instability point. Meanwhile the boson radius has its maximum value 
(27r) -1,/2 for the isotropic antiferromagnetic chains and decreases by changing in A and 
goes to zero at the ferromagnetic instability point. 

The maximum value of the spin excitation velocity occurs where K ~ 1/2 or 
A = — 1. It decreases monotonically by increasing K and reaches to the zero value at 
K = oo (A = 1). 
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The continuum limit of the Hamiltonian (|3j) is obtained by writing spin operators 
in terms of bosonic fields: 



\ —d x <f>+ ( - sm[V^K4>(x)}, 

V 7T 7T 



r+ = ^ e - i v / f^')((-l)« bsm[V^K<i>{x)} +c), 

(19) 

where a, 6 and c are non-universal real constant in the order of unity and depends on the 
parameter A [31]. In above transformations we have made the rotation — > (— l) n r^ 
and — > on the standard bosonic version of the spin operators which is found 
extensively in the literature. Using ( Jl9[ ), the effective Hamiltonian density is given as: 

n = v -[(dj) 2 + (d x( j>) 2 ] 

+ ^— sm[V^K<f)] - h e J f \[?-d x( j). (20) 
tioq V vr 

The mapped model is nothing but the massive sin-Gordon model with an additional 

topological term. Let us first consider the sin-Gordon model without the topological 

term h^f = or at the magnetic field value h = J af — if - As it has been discussed in 

many places, in the interval 1 < K < 2 the spectrum of the SG model contains soliton 

and anti-soliton with mass M.. The exact relation between the soliton mass and the 

bare mass h e ^ is as follows (35] 

M =C{K)h\ fS ^ K , (21) 



2-K 



0rT(| + i) 1 2T( 2 , 

By substituting K = 3/2 and h\ = SJaf i n (21) the soliton mass or equivalently the 
excitation gap is given as 

gapped SG model. This term creates a shift on the field 0. Actually, in the absence 
of the topological term the excitation spectrum of the model is gapped and the field 
4> is sticken on one of the minima where (sinv^Svr) = — 1. This one is correspond to a 
staggered AF order in the effective spin chain. Tuning the uniform field h e Q the number 
of particle is decreased and the minima is shifted. Challenging between two uniform 
and staggered field causes the system experiences a different phase when |/iq | > M. 
[36] . This phase transition is occurred at h~ where the excitation gap reopened. 
Consequently the width of the magnetization plateau is obtained as 

ht - K -f (* Jaf) 2 , (22) 
where h£ and h~ are the middle critical field and are given 

ht = J AF (l + j±—). (23) 
4 a 

Summarizing, implementing two continuum filed theories we found that the model 
has four critical fields. 



Space modulation on ID alternating spin- 1/2 model 



10 




Figure 3. (Color online.) Difference between the two lowest energy levels of the 
original Hamiltonian versus the magnetic field h, for chains with different lengths 
N = 16, 20 and exchange parameters Jp = 1.0, J af — 9/2 and S = 1/9. 

4. Numerical results 

To explore the nature of the spectrum and the quantum phase transition, we have used 
Lanczos method to diagonalize numerically chains with length up to iV = 28. 

First, we have computed the three lowest energy eigenvalues of chains with Jp = 1.0, 
different values of the length and antiferromagnetic exchanges. To get the energies of 
the few lowest eigenstates we consider chains with periodic boundary conditions. 

In Fig|3j we present results of these calculations for the exchanges Jp = 1.0, J af = 
9/2, 5 = 1/9 and chain sizes iV = 16, 20. We define the excitation gap as a gap on the first 
excited state. As it is seen in Fig. [3j in strong limit of antiferromagnetic exchange, this 
difference is characterized by the indistinguishable (within the used numerical accuracy) 
dependence on the chain length and shows an universal linear decrease with increasing 
magnetic field. At h = the spectrum of the model is gapped. Turning the magnetic 
field the energy gap decreases linearly with h and vanishes at h Cl . This is the first 
level crossing between the ground-state energy and the first excited state energy. The 
spectrum remains gapless for h ci < h < h~ , whereas the gap is reopened when h > h~ . 
After an increasing and a decreasing the spin gap goes to zero and vanishes at h^. 
More increasing in the field h > h^, the spectrum remains gapless up to the critical 
saturation field h C2 . Finally, at h > h C2 the gap is reopened and for a sufficiently large 
field becomes proportional to h. Oscillations of the energy gap in regions h Cl < h < h~ 
and h C2 < h < h£ are the result of level crossings in finite size systems. To find the 
critical fields we have used the phenomenological renormalization group (PRG) method 
[12J. The critical field values are given as follows: 

h ci =3.80 ±0.01, h C2 = 4.71 ± 0.01, 
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h c - = 3.99 ± 0.01, h c + = 4.58 ± 0.01. (24) 

To study the magnetic order of the ground state of the system, we start with the 
magnetization process. First, we have implemented the Lanczos algorithm on finite 
chains to calculate the lowest eigenstate. The magnetization along the field axis is 
defined as 

i N 

M z = j^J2( gs \ s j\ gs ^ ( 25 ) 

3=1 

where the notation (Gs\...\Gs) represents the ground state expectation value. In Fig. |1J 
we have plotted M z as a function of the magnetic field h, and for a chain with exchange 
parameters J F = 1.0, J A f = 9/2, 5=1/9 and different lengths N = 20, 24, 28 . As it is 
clearly seen in Fig. [4] besides the standard singlet and saturation plateaus at h < h Cl and 
h > h C2 respectively, we observe a plateau at M = ^M sat . Observed oscillations of the 
magnetization at h Cl < h < h~ and < h < h C2 result from the level crossing between 
the ground and the first excited states of this model in the gapless phases. To check 
that the mid-plateau is not a finite size effect, we performed the size scaling [37] of its 
width and found that the size of the plateau interpolates to finite value when iV — > oo. 
In the inset of Fig. [4jwe have plotted the magnetization on site, M? = (GS\Sj\GS), 
as a function of the site number "j" for a value of the magnetic field corresponding to 
the plateau at M z = 0.5M* at . To obtain an accurate estimate of the function Mf, we 
have calculated it for system sizes of iV = 12, 16, 20, 24, 28. The thermodynamic limit 
(N —> oo) of the finite size results are obtained by extrapolation method and used for 
plotting. As we observe the system shows a well pronounced modulation of the on site 
magnetization, where magnetization on odd bonds is larger than on even bonds. This 
distribution remains almost unchanged within the plateau for h~ < h < h c+ . 

By analyzing the numerical results on the energy gap (Fig. [3]), we found that the 
spectrum is gapped in absence of the magnetic field which is one of the properties of the 
Haldane phase. The Haldane phase can be recognized from studying the string order 
parameter. The string correlation function in a chain of length N, defined only for odd 
I, is EH 



2j+m 

Ostr(l,N) = -(exp{in £ S*}). (26) 

2j+l 

In particular, we calculated the string correlation function for different finite chain 
lengths. Since the present model has a SU (2) symmetry in the absence of a magnetic 
field, we only consider the z component of the string correlation function. In Fig. (|5|, 
we have plotted Ostr{l,N) as a function of h for the chain with exchanges Jp = 1.0 , 
Jaf = 9/2, 5 = 1/9 and lengths = 20,24,28. As can be seen from this figure, at 
h < h ci , the string correlation function Ostr{l, N) is saturated and the tetrameric chain 
system is in the Haldane phase. The Haldane phase remains stable even in the presence 
of a magnetic field less than h ci . 
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Figure 4. (Color online.) The magnetization along the field, M z , as a function of the 
applied magnetic field h, for chains with exchanges Jf — 1.0 , Jaf — 9/2, 6 = 1/9 
and lengths TV = 20, 24, 28. In the inset, the magnetization on site as a function of 
the site number j is plotted for a value of the magnetic field in the region of the to 
mid-plateau. 




o I- , , 



N=20 
N=24 
N=28 



Figure 5. (Color online.) The string correlation function, Ostr{l, N), as a function of 
the applied magnetic field h, for chains with exchanges Jp = 1.0 , Jaf =9/2, S = 1/9 
and lengths N = 20, 24, 28. 
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Figure 6. (Color online.) The AF-bond dimerization order parameter as a function 
of the applied field h, for chains with exchanges Jf = 1-0 , Jaf = 9/2, 5 = 1/9 and 
lengths N = 20, 24, 28. 



An additional insight into the nature of different phases can be obtained by studying 
the correlation functions. We define the following weak and strong bond dimerization 
order parameters; 

4 N/2 

d r=N £ (GS\S 2 j-i ■ S 2 j\GS), (27) 

jodd=l,3,5,... 

and 

4 N/2 

d r = N £ (GS\S 2j ^.S 2j \GS), (28) 

Je.ve.ri — 2,4,6,... 

where summations are taken over the weak and strong antiferromagnetic-bonds. In 
Fig. ^ we have plotted d™ and d s r versus magnetic field h for chain of lengths 
N = 20,24,28 with the exchange parameters Jp = 1.0, Jaf = 9/2 and 5 = 1/9. 
As it is seen from this figure, at h < h Cl spins on all antiferromagnetic-bonds are in 
a singlet state d™ = d s r ~ —0.75, while at h > h C2 , d r is equal to the saturation value 
d™ = d s r ~ 1/4 and the ferromagnetic long-range order along the field axis is present. 
However, in the considered case of strong antiferromagnetic exchanges (J^ F S> Jf) and 
high critical fields quantum fluctuations are substantially suppressed and calculated 
averages of on-antiferromagnetic-bond spin correlations are very close to their nominal 
values. 

On the other hand, for intermediate values of the magnetic field, at h Cl < h < h C2 
the data presented in Fig. ^ gives us a possibility to trace the mechanism of singlet-pair 
melting with increasing magnetic field. As it follows from Fig. ^ at h slightly above h Cl 
spin singlets pairs start to melt in all antiferromagnetic bonds simultaneously and almost 
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with the same intensity. With further increase of h melting of weak antiferromagnetic- 
bonds gets more intensive, however at h — h~ the process of melting stops. As it is 
seen in Fig. ^ weak antiferromagnetic-bonds are polarized, however their polarization 
is far from the saturation value d™ ~ 0.15, while the strong antiferromagnetic-bonds 
still manifest strong on-site singlet features with d s r ~ —0.65. Moratorium on melting 
stops at h = however for h > h+ strong antiferromagnetic-bonds start to melt 
more intensively while the polarization of weak antiferromagnetic-bonds increases slowly. 
Finally at h — h C2 both subsystems of antiferromagnetic-bonds achieve an identical, 
almost fully polarized state. Note, that the almost symmetric fluctuations in on- 
antiferromagnetic-bonds correlations, increase in d™ at h < h~ decrease in d s r at h > h£ 
reflect the enhanced role of quantum fluctuations in vicinity of quantum critical points. 

In our previous paper [T3], we introduced a mean field order parameter which can 
distinguish a gapless LL phase from the other gapped phases. This order parameter is 
the F-dimer order parameter which is defined as 

P F = Re(S^ n St n+1 ). (29) 

The F-dimer order parameter has a considerable value in the Liittinger liquid phase 
and behaves differently in the other gapped phases. The effects of a small value of 
space modulation on this parameter has been shown in Fig. ([7]). We have plotted Pp 
versus magnetic field h for chain of lengths N = 12, 16, 20, and 24 with the exchange 
parameters Jp = 1.0, Jaf = 9/2, and 5 = 1/9. As it is seen, in the Haldane phase, 
h < h Cl , the quantum fluctuations suppress the ferromagnetic correlations and the F- 
dimer parameter is close to zero-value. Right after the first critical field, the F-dimer 
parameter increases rapidly up to the second critical field. In the intermediate region, 
h c - < h < h c +, the F-dimer parameter shows a non-zero plateau which behavior is the 
same as the other parameters. By more increasing field, the F-dimer decreases and goes 
to zero at the saturation critical field h = h C2 . The intermediate region is a gapped 
phase and we expect the zero-value for the defined LL parameter, Pp. However the gap 
does not affect on the behavior of Pp and a plateau appears in the curve. 



5. Pair- wise entanglement 

In this section we focus on the entanglement of two spins in different phases of the 
system. The entanglement which has no a classical counterpart is employed to study 
the quantum correlations of different states. Concurrence is a measure of the bipartite 
entanglement which is defined as following [351 HO] 

C lm = 2max{0,C i l 1 J,Cl 2 J}, (30) 

where 

/] M z — M z 

1 2 m ) 2 , 
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Figure 7. (Color online.) The F-dimer order parameter as a function of the applied 
field h, for the chains with exchanges Jf = 1-0 , Jaf = 9/2, and S — 1/9 and lengths 
N = 12, 16, 20, and 24. The appearance of a non-zero plateau in the curve is clear. 

Cf = yj(9% + 9Z) 2 + W " O 2 

- V^ + ^) 2 -( j 3 T (31) 

and = (S^S^) is the correlation function between spins / and m. The numerical 
Lanczos results on the concurrence for the ID tetrameric spin-1/2 model have been 
shown in Fig. (J8|. We have plotted the entanglement of two spins which are located 
at the same strong, weak, and ferromagnetic bond versus h, chain length N = 28 and 
the exchange parameters Jp = 1.0, Jaf = 9/2, and 5 = 1/9. In the Haldane phase, 
h < h Cl , the spins on all antiferromagnetic bonds make a singlet state. In this state 
which is a maximally entangled state, gf£ = gf^ = gf^ = C s = C w = 1, and C F is 
zero. For h > h Cl the values of C s and C w fall down with increasing of the magnetic 
field. Indeed the quantum correlations of the two spins with strong antiferromagnetic 
(SSA) and weak antiferromagnetic (SWA) interactions are decreased by increasing the 
magnetic field. However an enhancement on the entanglement of the two spins with 
ferromagnetic interaction (SF) is observed. It means that the magnetic field increases 
the quantum correlations of the two spins which are interacting ferromagnetically. This 
is a dual effect of the magnetic field in which increases the quantum correlations of two 
SF and decreases the quantum correlations of two SWA and two SSA. In the intermediate 
gapless region h Cl < h < h~ , the quantum correlations of SWA and SSA diminish down 
to h — h~ and the concurrences C w and C s reduce to ~ 0.1 and ~ 0.9 respectively. 
However, the quantum correlations of SF grows up to the critical field h c _ and the 
entanglement reach to the value ~ 0.25. At the plateau state, the gap of the system 
is re-opened and a plateau emerges in the curve of concurrences. In the intermediate 
gapped phase the values of the concurrences C s , C w , and C F are 0.9, 0.1, and 0.25, 
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Figure 8. (Color online.) (a)-The concurrence as a function of the applied field h, for 
chains with exchanges Jp = 1.0 , J^f = 9/2, and 5=1/9 and length N = 28. Inset 
plot: The concurrence between two spins on strong bond as a function of h for different 
chain sizes N — 16, 20, 24, and 28. (b)- Concurrence of two spins versus separation 
distance m. 



N =28 




respectively. Indeed in the plateau state there are three types of quantum correlations 
in the system. These correlators are the source of the mid-plateaus in the different 
parameters of the system such as magnetization, bond dimer and F-dimer parameters. 
Indeed all of these quantum correlators are exist only at the mid-plateau states. In the 
full saturated state all of them disappear and the entanglement of the state is exactly 
zero. 

To see the finite size effects, we have plotted in the inset plot of Fig. |8](a), the 
concurrence between two SSA as a function of h for different chain sizes N = 16, 20, 
24, and 28. It can be seen that there is not size effect on the numerical results and the 
concurrence behaves as thermodynamic limit in the gapped regions. 
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To get more intuition on the mid-plateau state, we have computed numerically the 
entanglement between the spins with different separation distances. The concurrence 
between a SSA and an arbitrary spin, say S m has been plotted versus the separation 
distance m, in Fig. [8]^b). The entanglement of two such spins is decreased by increasing 
m and goes to zero at the finite value of m = 7. The same behavior is also observed for 
the entanglement between a SWA and the spin S m with respect to m. The entanglement 
between SWA and S m goes to zero at m — 9. It means that, in the plateau state the 
range of the quantum correlations between a SWA and S m is longer than the range of 
a SSA and S m . This is the other feature of the mid-plateau state. 

It is also remarkable that from our numerical results we found that the concurrence 
of two spins that are not on the same bond is equal to zero in the gapped Haldane 
and saturated ferromagnetic phases which is in complete agreement with the analytical 
results. 



6. conclusion 



In this paper, we have focused on the magnetic properties of a bond-alternating 
antiferromagnetic-ferromagnetic spin- 1/2 Heisenberg chain. Using two analytical 
approaches and a numerical method we have studied the effects of an external magnetic 
field and a space modulation on the ground state properties of the system. In the limit 
where the AF couplings are dominant, we mapped the model (fil) to an effective XXZ 
Heisenberg chain with anisotropy parameter A = 1/2 in the presence of effective uniform 
and staggered longitudinal magnetic fields. This model has different quantum phases 
which are distinguished by four critical fields. To find the critical fields we employed two 
field theoretical approaches such as nonlinear sigma model and bosonization, and the 
numerical exact daigonalization Lanczos method. Working on the spin coherent states 
representation, we mapped the model Q to a nonlinear sigma model with an additional 
topological term. The topological term is dependent on the space modulation 5 and the 
parameter a = j^. For any values of a and 5, the Hamiltonian mapped to a non- 
integrable 0(3) NLcrM with in the interval [0, 7r]. In the absence of a magnetic field 
the model is always gapped and the gap value depends to a and 5. By analyzing our 
NLctM in presence of the magnetic field, we obtained only the two critical fields. To 
dominant this vacancy and to find the other two critical fields we also bosonized the 
effective Hamiltonian of XXZ chain ([3]). Our bosonization procedure showed that the 
width of the mid-plateau is a function of 5Jaf- It has been shown that the plateau- 
width, scales as a power low with exponent value 2 and vanishes at 5 = 0. 

Moreover, we implemented the Lanczos method to numerically diagonalize chains 
with finite length up to iV = 28. Using the exact diagonalization technique, we 
calculated the energy gap, magnetization, the string, the F-dimer, and the bond dimer 
order parameters and various correlation functions for different values of the external 
magnetic field. In good qualitative agreement with our analytical results, we showed 
clearly that a space modulation on the antiferromagnetic exchanges leads to generation 
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of a gap in the excitation spectrum of the system and correspondingly a magnetization 
mid-plateau at = \. We found that a non-zero plateau also creates in the plot of 
F-dimer, bond dimer order parameters. 

To get more physical insight on the mid-plateau state we also investigated the 
pair-wise entanglement between two different spins of the system. As a measure of 
entanglement the concurrence between two arbitrary spins computed as a function of 
magnetic field. A plateau is also appeared in the concurrence at the middle gapped 
state. In the plateau state there are three types of quantum correlations in the system, 
in the plateau state the range of the quantum correlations between a spin on the weak 
antiferromagnetic bond and a S m is longer than the range of the quantum correlations 
of a spin on the strong antiferromagnetic bond and S m . 
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